Future Probes of the Primordial Scalar and Tensor Perturbation Spectra: 
Prospects from the CMB, Cosmic Shear and High- Volume Redshift Surveys 
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Detailed study of the scalar and tensor perturbation spectra can provide much information about 
the primordial fluctuation-generator, be it inflation or something else. The tensor perturbation 
spectrum may be observable through its influence on CMB polarization, but only if the tensor-to- 
scalar ratio, r = T/S, is greater than about 10^^. The tensor tilt can be measured with an error of 
a{nT) that decreases with r from 0.1 at r = 0.001 to 0.02 at r = 0.1. Current CMB constraints on 
the scalar perturbation spectrum can be improved by higher-resolution CMB observations and/or 
by tomographic cosmic shear observations. These can both shrink errors on the tilt (ns) and 
running (n'g = dns / dink) to the 10"'^ level. Stunning as these results would be, it may become 
very desirable to improve upon them an order of magnitude further in order to study the expected 
departures from n'^ — 0. Such improvements are likely to require observation of three-dimensional 
clustering over very large volumes. Unfortunately, to get down to the 10~* level will require a sparse 
spectroscopic redshift survey with about 10^ galaxies spread over a volume less than but comparable 
to that of the observable Universe. 



INTRODUCTION 



Inflation is doing remarkably well with respect to ob- 
servations. The evidence that structure in the Universe 
formed from adiabatic, nearly scale- invariant fluctuations 
is very strong 0]. The mean spatial curvature has been 
determined with high accuracy and is consistent with 
zero The perturbations are highly Gaussian 

These are all predictions of inflation, of varying degrees 
of robustness. Topological defect models, formerly in- 
flation's chief competition, have been strongly ruled out 
(e.g., i). 

Despite this great success we still have little under- 
standing of the physics that led to these initial perturba- 
tion spectra. There are many different models of infla- 
tion. As Michael Turner said in his talk at this meeting, 
we have a Landau-Ginzburg theory, but we are missing 
the underlying BCS theory. 

To discriminate among inflationary models, to assist 
in the quest for a BCS theory of inflation, and perhaps 
to discriminate between inflation and alternatives (such 
as we heard from Steinhardt and from Durrer), we want 
to take the determination of the primordial tensor and 
scalar perturbation spectra to a qualitatively new level of 
precision. This point was also emphasized in Albrecht's 
talk. 

I will first briefly review the current constraints on the 
scalar perturbation spectra before turning to the future. 
The discussion of future experiments will start with the 
constraints on the tensor perturbation spectra from CMB 
experiments. We will then look at how current CMB 
data can be complimented cither with higher-resolution 
CMB data and/or tomographic cosmic shear. Finally, I 
consider results that can in principle be achieved by very 
high- volume spectroscopic redshift surveys. 



II. NON-ZERO RUNNING? 

A question of great relevance to inflationary models 
is the significance of evidence for running of the spec- 
tral index. The evidence from CMB data alone is very 
weak: n'g = dus/dlnk = -0.055 ± 0.038 0]. Combining 
CMB data with the Croft et al. matter power spectrum 
inferred from high-resolution observations of the Lya for- 
est results in n'g = -0.031 ± 0.017 j^]. Other authors H 
working with the same datasets have since found much 
looser constraints on n'g. The looser bound on n'g is due 
to their marginalization over the mean ionizing flux as 
a function of redshift. The importance of marginalizing 
over this parameter, which leads to a large degeneracy 
between spectral index and amplitude, was pointed out 
in 0. 

We had much discussion at this meeting about con- 
straints from combining CMB data with low-resolution 
spectra from the Sloan Digital Sky Survey (SDSS), as 
presented by Hui and by Seljak. Although lower reso- 
lution, there are thousands of these spectra allowing for 
very small statistical errors. Possible systematic errors 
have not yet been understood and controlled well enough 
to allow for any firm detections of n'g ^ 0. 

To summmarize the current situation: there is no 
strong evidence for r\!g ^ 0, but very interesting results 
may be coming soon from SDSS quasar spectra combined 
with WMAP. 



III. CMB AND TENSOR PERTURBATIONS 

One can decompose a polarization pattern on the 
sky into curl-free modes (E modes) and divergence-free 
modes (B modes) ^, 8]. Each scalar three-dimensional 
Fourier mode only has one direction, that given by the 
wavevector k, and therefore only leads to E modes. Ten- 
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sor Fourier modes (gravitational waves) have the direc- 
tion given by k and also an orientation given by the po- 
larization of the gravitational wave. They can therefore 
generate B modes in addition to E modes. Since, at least 
in linear perturbation theory, scalar perturbations do not 
generate B modes, whereas tensor perturbations do, the 
B mode has been proposed as a means to detect the ten- 
sor perturbations. 

At second order in perturbation theory we cannot solve 
for just one Fourier mode at a time and then sum up the 
resulting solutions; the evolution of one Fourier mode 
is affected by the presence of others. This mode-mode 
coupling leads to the generation of B modes even from 
scalar perturbations The dominant second order 

effect is gravitational lensing. The lensing-induced B 
modes can obscure the gravitational wave contribution 
to the B modes. If the lensing-induced B modes are not 
cleaned from the B-mode map, one can only detect the 
tensor signal (at 3ct) if r > nim = 1-7 x 10""*, 7 x 10"^ 
or, 2.4 X 10^^ for r — 0.05,0.1 or 0.2 respectively. 

Gravitational lensing leads to off-diagonal correlations 
in Fourier space 10, 11, 12, 13]. These can be used to re- 
construct the lensing potential ||l3,Eli ■ With the lensing 
potential thus reconstructed, the maps can be unlensed. 
A perfectly unlensed map would have no scalar B modes 
in it. Unavoidable imperfections in the lensing potential 
reconstruction mean the unlensed map will have some 
residual B mode, even in the absence of tensor pertur- 
bations. These residual B modes prevent detection of 
gravity waves unless r > riim = 1.4 x 10~^,6 x 10^^ or 
2 X 10"^ for r = 0.05,0.1 or 0.2. These limit calculations 
were first done by [1^1 and p^ . 

As shown in [33, if r > 0.01 then we can learn some- 
thing about the shape of the tensor B mode power spec- 
trum. The possible constraints on ut as a function of r 
are shown in Fig. 1. 

Tensor spectra from inflationary models with a single 
field in slow roll obey the 'consistency equation', r — 
—briT 18]. Note that, more generally, r = f(yL\)nT and 
that /(O) ~ -7 and /(0.7) ~ -5 From the figure 

we can see that for r > 0.03 we can make a significant 
test of the consistency relation. The consistency equation 
applies for all single-field models to first order in the slow 
roll parameters. 



IV. CMB AND SCALAR PERTURBATIONS 

Constraints to the scalar perturbation spectra can be 
improved by pushing to higher angular resolution than 
the ~ 13' of the WMAP's highest frequency channel. It 
can not be improved arbitraily though because of the ex- 
ponential suppression of power that sets in at the Silk 
damping scale, due to photon diffusion during recombi- 
nation. High sensitivity can to some degree fight against 
this exponential cutoff, but not to indefinitely high £. 

Here we show results expected for three different ex- 
periments (with parameters specified in Table I for the 




FIG. 1: Error on nr + r/4.8 from a full-sky, no-noise ex- 
periment (solid line) and for varying angular resolutions (top 
panel) and varying noise levels (bottom panel), as a function 
of the tensor-to-scalar ratio, r. The inflationary consistency 
equation is ut + r/4.8 = 0. Values of r > 0.71 are ruled out 
by temperature power spectrum data. From Song & Knox, in 
preparation. 

Experiment ij^^x ^mkx (GHz) 6t At Ap 



Planck 2000 3000 100 9.2' 5.5 00 

143 7.1' 6 11 

217 5.0' 13 27 

SPTpol (/sky = 0.1) 2000 3000 217 0.9' 12 17 

CMBpol 2000 3000 217 3.0' 1 1.4 



TABLE I: Experimental specifications. 

amplitude of the primordial gravitational potential power 
spectrum, P^{kf), ns, and n'g where 

In [Plik)] - In [Pl{kf)] + {ns{kf) - 3) In (k/kf)+n's[\n {k/kf)f. 

(1) 

We also include the expected constraints on w and m^. 
The errors in Table II are what would result from a si- 
multaneous fit to these parameters plus flbh^, flmh^, 
(the angular size of the sound horizon), the Helium mass 
fraction and r. We assume the mean spatial curvature is 
zero. See for more details. 

Are these highly precise measurements of ns and n'g 
valuable? The difference ng — 1 is first order in the slow- 
roll parameters and n'g is second order. Thus we expect 
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Experiment 


(cV) 


VJx 


In Pi 


ns 


n's 


T-'l ct Yi r'lr 

SPTpol 
SPTpol + Planck 
CMBpol 


0.14 
0.11 
0.082 
0.031 


0.28 
0.34 
0.22 
0.088 


0.016 
0.018 
0.016 
0.011 


0.0074 

0.01 
0.0057 
0.0024 


0.0032 
0.0057 
0.0027 
0.0014 







TABLE II: Standard deviations expected from Planck, SPT- 
pol and CMBpol. From 21]. 



\n'g \ to be on the order of (715— 1)^. If rig— 1 — 0.03, which 
is perfectly consistent with present data, then we would 
expect to see non-zero n'g at about the 10"'^ level. It 
would be tremendously exciting to actually confirm this 
expectation. Note though that there are models with 
\n'g\ ~ 715 — 1 22] so this is not a firm test. Nevertheless, 



finding ij{n'g) >> ns 
with inflation 1231. 



1 would be difficult to reconcile 



V. CMB+COSMIC SHEAR 

Adding cosmic shear to CMB data can improve con- 
straints of cosmological parameters, as was studied in the 
first CMB+cosmic shear forecasting paper j23]. A num- 
ber of other studies have followed, extendingto use of 
photometric rcdshifts to allow tomography |25l ] and to 
examine constraints possible on Wx in addition to the 
neutrino mass 0, [23, 113 . Here we focus on errors on 
ns and n^, though also make some comments about con- 
straints possible on and Wx- 

The shapes of galaxies are distorted by gravitational 
lensing. If we knew the shapes of the galaxies in the 
absence of lensing we could then infer the lensing conver- 
gence, K, from the lens-induced alteration of the shape. 
Although we have no way of knowing the unlensed shape 
of an individual galaxy, we expect that an average over a 
large number of galaxy images would be perfectly cylin- 
drical. The average departure from a completely cylindri- 
cal shape can thus be used as a measure of the amount 
of lensing. From proper analysis of the average galaxy 
shape in some pixel, the lensing convergence, k, can be 
recovered. 

We model the data n (the convergence vector with in- 
dex running over m and redshift bin 2;) as a Gaussian 
random field with zero mean and contributions to the 
total variance from both signal and noise: 



(kk^) = S + N. 
The signal covariance has the structure 



Sir, 



(2) 



(3) 



and the noise covariance is entirely diagonal: 

Nlmz.J'm'z' = NiSii>6mm'5zz' ■ 



(4) 



The average ellipticity of galaxies, jrms, leads to an 
error in the estimated real-space convergence map with 
variance 7rms/^z"^ where N^^^ is the mean number of 
galaxies in each pixel in redshift bin z. This error is 
uncorrelated from pixel to pixel and from redshift bin 
to redshift bin. It leads to an error in the Fourier- 
transformed convergence with variance |29| 



(5) 



For specificity we take ilpix = 11 sq. arcmin and 
7rms = 0.2. The mean number of galaxies in each pixel 
we take to be 



14.2zi-iexp[-(z/1.2)i-2] 



(6) 



for Az = 0.4. This is a fit to the extrapolation from 
Nagashima et al. 2002 assuming a limiting magnitude in 
R of 26 (Tony Tyson, private communication). We also 
assume that half of these galaxies in the 1.2 < z < 2.5 
range we will not be able to use because we will not be 
able to get sufficiently accurate photometric redshifts. 

As the pixel size increases, the number of galaxies in- 
creases and the statistical error drops (equivalently, in 
Fourier space, PNi increases with I since Ni is constant). 
We do not expect any observations to keep the system- 
atic errors smaller than 10~^ on scales above 10'. With 
our modeling of the statistical errors, the rms statisti- 
cal error on the shear drops below 10"'* on scales of 5 
degrees. For this reason, we quote our results without 
/ < Ic where Ic = 180/5 = 36. 

The auto- and cross- angular power spectra for the 
convergence, Cf^ , are given by (e.g., 0) 
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drrWir)^ (r)A|(fc,r) 



(7) 



with / — kr and 

1 2 ,(r'-r) 



W'{r) = / dr'^ —!-N,{z')dz'/dr'Q{r' ~ r) 

^ Jr(z,,) r' 

(8) 

where zn — Zi — A.z/2, Zi2 ~ Zi + Az/2 and n.i is 



r{zn) 



dr' Ni{z')dz' / dr' 



(9) 



and we have assumed that the galaxies are evenly dis- 
tributed through the redshift bin. 

In addition to the convergence, k, there are also the 
lensing potential, and the deflection angle d. They are 
related in real space by d = and k = V'^0/2. Their 
power spectra are related by 

2(^K(.)«(.')/^ = ;(Z+l)Cf '^)''(^V(27r) = Z2(;+i)2cf (^)^(^'V(27r) 

(10) 
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In Fig. 2 we see the convergence power spectra for four 
redshift slices: z = 0.6, z = 1, z = 2.2 and z = 1100. Er- 
rors on z = 1100 are those expected from Planck. Errors 
on the others are those expected from the 30,000 sq. de- 
gree LSST survey described above with bins in redshift 
of width Az = 0.4. The convergence rises monotoni- 
cally with redshift since the lensing contribution from 
each redshift interval adds incoherently to that of the 
previous redshift interval. 



io-« 



10- 



5^ io-» 



io-» 



10-10 



FIG. 2: Convergence power spectra for four redshift slices: 
z = 0.6, z = 1, 2 = 2.2 and z = 1100. Errors on z = 1100 are 
those expected from Planck. Errors on the others are those 
expected from a 30,000 sq. degree survey to mjj = 26 with 
bins in redshift of width Az = 0.4. 




In Fig. 



Cf 



: ^z' z' 



for 



plot the correlation, r 
z = 0.6, 1 and 2.2 with z 



zz 

I 



fixed 



to 1100. The correlation increases with increasing red- 
shift as the window functions become better matched. 
In other words, the correlation rises with increasing red- 
shift since as z goes towards 1100 a greater fraction of 
the lensing structures is common to both. 

By differencing the auto power spectra with appropri- 
ate weightings one can isolate how much lensing power 
is coming from each redshift range, and therefore the 
power spectrum of the matter in each redshift range. To 
forecast errors on the parameters that govern the matter 
power spectrum, we need not do this differencing. Our 
statistical model of the data is completely specified by 
the auto and cross-power spectra and we need only un- 
derstand our errors on these and how they vary as we 
vary the model parameters. 

To forecast errors we calculate the Fisher matrix 
1. 




1000 



for 



FIG. 3: The correlation, r,"' = C^' / ^/cf 
0.6, 1 and 2.2 with z' fixed to 1100. Errors are those expected 
from Planck plus a 30,000 sq. degree survey to niR — 26 with 
bins in redshift of width Az — 0.4. 



I,zl,z2,z3,z4: 

where the subscript , p means differentiation with respect 
to parameter Up and 



W = (S + N)-\ 



(13) 



The inversion of W is numerically tractable because of 
its block-diagonal structure: 



mz,l' m' z' 



= 5m.m'5z 



(14) 



J~ Ti' 



-Tr [S,pWS,p-W] 



(11) 



Since photometric redshifts are completely unproven 
for z > 3 we make no use of galaxies at z > 3. To 
reduce sensitivity to the non-linear evolution which we 
may not have modeled well enough we restrict ourselves 
to / < 1000. Results are shown in Table III. 

There is not much to be gained by going deeper, at 
least with our z-independent I cutoff at I = 1000. The 
reason is that the number density of galaxies is high 
enough in each redshift bin that the Ci measurements 
are sample-variance limited for / < 1000. For the lower 
redshift bins there is not much to be gained by going to 
higher I because that puts us into the fully non-linear 
regime which has little to no shape information. Also in 
the non-linear regime our formalism for forecasting the 
uncertainty breaks down. But for higher redshift bins 
one can in principle push to higher I and gain by it. The 
I that divides linear from non-linear is higher at high red- 
shift for two reasons: l)a given k projects into a higher 
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Experiment 


m„ (eV) 






ns 


n'g 


WMAP 
Planck (no lensing) 
WMAP+LSST 
Planck+LSST 


1.2 
0.32 
0.085 
0.041 


2.4 
0.45 
0.023 
0.019 


0.061 
0.016 
0.023 
0.012 


0.055 
0.0077 
0.0094 
0.0042 


0.018 
0.0035 
0.0027 
0.0021 







TABLE III: Standard deviations expected from WMAP (after 
2 years), Planck (without reconsruction of Cf''), WMAP + 
LSST and Planck (with reconstruction of Cf'*) + LSST. 



I, and 2) the linear to non-linear transition is at higher 
k. In principle our / cutoff should be z-dependent but we 
use an /-independent one for simplicity. 

There is also not much to be gained by going to higher 
redshift. There are two reasons for this. First, the an- 
gular diameter distance does not change very much over 
this range of redshift so we are not probing significantly 
different spatial scales. And second, we are already close 
to the sample variance limit so the noise reduction due 
to the extra galaxies is not useful. Again, this conclusion 
could change if we allowed a z-dependent I cutoff that 
followed the linear/non-linear transition. 

In addition to the precision determination of the pri- 
mordial perturbation spectrum, we can also measure the 
dark energy equation-of-state parameter Wx and the neu- 
trino mass to extraordinary precision. With data that 
can measure Wx to 0.02 we could clearly also place strong 
constraints on its evolution. Note, however, that rrii, is 
no better determined than it can be by the CMB alone 
(see Table II) with CMBpol. 

In the first cosmic shear + CMB forecasting papers, 
p^ l and [2^, it was shown that combining Planck and 
all-sky cosmic shear could achieve neutrino mass errors 
of 0.04 eV without tomography and 0.02 eV with to- 
mography, with Wx fixed to —1. We find similar results 
when we fix Wx to -1. More recently [23|, [1^ and 
also studied constraints on neutrino mass and/or Wx that 
can come from tomographic cosmic shear surveys. These 
studies are done with all other cosmological parameters 
held fixed. The justification is that these parameters will 
be determined by CMB observations. We find, however, 
that the remaining uncertainties in cosmological parame- 
ters, even given the expected CMB data, are significant. 
This remaining uncertainty greatly increases errors on 
in particular, and to a lesser degree on Wx- 

It may not be possible to control systematic errors at 
the requisite levels on large angular scales. We coarsely 
simulate the effect of systematic errors on large angular 
scales by imposing a low I cutoff at I = 36; i.e., we throw 
out all cosmic shear information at Z < 36 (though we 
keep the CMB-determined convergence map all the way 
down to Z = 2). The resulting error increase from this 
cut is small. Most of the weight is coming from higher 



I. Although we need full sky to get these exquisitely 
small errors, we only need full sky to beat down sample 
variance, not to accurately determine flucutation power 
on large angular scales. 



VI. CMB + HIGH- VOLUME REDSHIFT 
SURVEYS 

The errors on n'g from WMAP-I-LSST or from higher- 
resolution CMB observations, are just at the level of 
{ns — 1)^ if ns — 0.95. These experiments may very 
well leave us only with tantalizing one or two a determi- 
nations of n'g. We may want to do better. 

Better CMB experiments will not do it. Better tomo- 
graphic cosmic shear experiments probably will not do 
it. We can understand why with a simple expression for 
the error on a power-spectral index, n, given errors on 
the power Pi and P2 at two different scales, fci and /C2: 



a^{n) 



<T^(lnPi)+g2(lnP2) 



(15) 



First note that gaining dynamic range (increasing fc2/^i) 
only reduces the error on n logarithmically. Extend- 
ing dynamic range is not a way to dramatically im- 
prove the constraints on n. The only way to dramati- 
cally reduce a(n) is to make better measurements of the 
power over around a decade (or more) in wavenumber. 
2-dimensional surveys can not do this. Even with 9 un- 
correlated measures of fluctuation power in redshift bins 
between zero and last-scattering we can reduce the CMB 
alone errors by at most a factor of 3. 

Thus we explore the constraining power of three- 
dimensional measurements of clustering. In particular, 
we consider a spectroscopic redshift survey. See [sJl for 
more discussion of future redshift surveys. For simplic- 
ity, we ignore effects of redshift distortions and scale- 
dependent galaxy bias, which should be included in a 
more detailed study. Note that very deep cosmic shear 
surveys may be invaluable for calibrating the galaxy- 
density relationship for these spectroscopic surveys (in 
order to constrain the bias), and for providing the tar- 
gets for spectroscopic follow-up. 

For a survey of volume V , the error on the power in a 
band of width Ak is (3^ i3^] : 



(16) 



where Nk is the number of independent samples of (5(k) 
with k - Afc/2 <\k\ < k + Ak/2 and is given by 



2.. V 



(27r)3 



k^Ak/{2TT^)V. (17) 



The mean number density of galaxies, n, determines the 
shot-noise contribution to the variance, l/(n&^) where b 
is the galaxy bias such that Pg{k) = b'^P{k). 
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An optimal survey with fixed number of galaxies Ng = 
nV (optimal from a solely statistical error point of view) 
will have a volume such that 1/n — b^P{k) on the scale 
of interest where b is the bias of the galaxies. In this case 



APjk) 

P{k) 



Ngb^A^ik) 



Ak 



-1-1/2 



where Ng — nV is the number of galaxies in the survey 
and A2(fc) = fc3p(fc)/(27r2). 
Solving for Ng we find 



Ng = 1.5 X lO-* 



10- 



n 2 r 



AP/P 



0.5 



62A2(fc) 



Ak 



(19) 



This is a lot of galaxies. And unfortunately, most of them 
are at high redshift (if we want AP/P = 10^"*). If the 
volume is spherical with coordinate distance radius of R 
then 



R = 7h-^Gpc 



0.2hMpc- 



(fc/Afc)i/3 



10- 



AP/P 



2/3 



(20) 

For reference, the distance to the horizon is about 14 Gpc 
[k,34]. 

One might conceivably do better by extending to 
higher k than the fiducial 0.2 h Mpc~^. The main ben- 
efit is the reduced volume one has to sample {V oc k~^, 
see above equation) to achieve the same value of AP/ P. 
The higher this is pushed, however, the more one must 
worry about scale-dependent bias. And when one wants 



AP/P as small as 10"'^ even a tiny amount of unknown 
scale-dependence can lead to highly significant system- 
atic error. 



VII. CONCLUSIONS 

If inflation happened at sufficiently high energy scales 
(1^1/4 > 2 X 10^^ GeV) and nature is kind to us with 
respect to astrophysical foregrounds (dust, synchroton 
radiation, etc.) then we can detect the influence of tensor 
perturbations in the CMB polarization. If the energy 
scale is slightly higher, we will be able to verify (or rule 
out) the inflation consistency equation. 

For measurement of the scalar perturbation spectrum, 
constraining ns and n'g to better than 10"'^ may be im- 
possible. This level can be reached with a post-Planck 
CMB polarization mission or by combining WMAP or 
Planck with all-sky tomographic cosmic shear observa- 
tions. Spectroscopic redshift surveys, even with benign 
assumptions about galaxy bias, are unlikely to improve 
our ultimate constraints on ns and n'g. 
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